The Ground State of a Two-dimensional Lattice System with a Long-range 
Interparticle Repulsion. Effective Lowering of Dimension. 

A. A. Slutskin, V. V. Slavin, and H. A. Kovtun. 

Institute for Low Temperature Physics and Engineering, 47 Lenin Ave., Kharkov, Ukraine. 

(*) 

It has been shown that effective lowering of dimension underlies ground-state space structure and 
properties of two-dimensional lattice systems with a long-range interparticle repulsion. On the basis 
of this fact a rigorous general procedure has been developed to describe the ground state of the 
systems. 



PACS number(s): 64.60. Cn, 64.60.-i 

The two-dimensional (2D) lattice systems with a long- 
range interparticle repulsion (LSLRIR) are of great in- 
terest as they have important nhysical applications. One 
of them is the adatom systemsa with a sufficiently strong 
interaction between the particles and the substrate. An- 
other field to apply the 2D LSLRIR model is a 'frozen' 
electron phase in 2D or layered narrow-band conduc- 
tors with a long-ranged electron-electron repulsion. This 
state results from a suppression of a narrow-band elec- 
trons' tunneling between host-lattice sites which is pro- 
duced by the mutual electron repulsiono and, therefore, 
differs principally from a Wigner crystal. (It arises if 
t < Su = (a/r)u, where t is the electron bandwidth, Su is 
the typical change in the energy of a narrow-band elec- 
tron as it hops between the host-lattice sites, a is the 
host-lattice spacing, r is the mean electron separation, u 
is the mean Coulomb energy per electron) . Layered con- 
ductors such as cuprates and polycrystal boundaries of 
the non-linear electroceramic materialscl as well as art 2D 
conductorsa appear to be most favorable for the electron 
'freezing'. 

At t <C Su the 2D - 'frozen' - electron phase ground- 
state space structure is much the same as that of the 
adatom systems. As far as we know, even in this limit 
only the one-dimensional (lEOpLSLRIR have been-,stud- 
ied adequately (ground stateETQ, thermodynamicsEl) . As 
to the 2D LSLRIR, neither its ground state nor thermo- 
dynamics have been explicated. Here we develop a rigor- 
ous procedure to describe, in the zero-bandwidth limit, 
the ground state of the 2D LSLRIR with an isotropic 
pair potential of the interparticle repulsion, v(r) (r is a 
distance between interacting particles). The basis of our 
consideration is zero-temperature effective lowering of di- 
mension in the 2D LSLRIR, which we have found to take 
place irrespective of a form of v(r) (within the limits of 
weak and physically reasonable restrictions)) and for any 
filling factor, p — N/AT (N and J\f are the total num- 
bers of the particles and host-lattice sites respectively, 
N, Af — > oo). We take the term to mean that the ground 
state of the 2D LSLRIR is an effective ID LSLRIR whose 
"particles" are strip-like periodic structures on the 2D 
host lattice. This enables to describe the ground gtate 
analytically, using results of the ID LSLRIR theorjoQ. 

In the limiting case under consideration the Hamilto- 



nian, H, of the 2D LSLRIR is of the form 

where r are radius vectors of sites of the host lattice; it is 
assumed to be a triangular one with a = 1; n(r) = or 1 
is the number of particles at a given site; the sum is taken 
over the whole host-lattice; function v(r) is assumed to 
decrease faster than r~ 2 as r — > oo and to be changed 
substantially over distances ~ r. 

It is reasonable to build up the 2D LSLRIR ground 
state description around the simplest 2D structures 
which are stable triangular 2D crystals with one particle 
per cell and p = p q = 1/q 2 (i.e. their primitive transla- 
tion vectors have a form qa, where a are those of the host 
lattice). We will call them "q-crystals" . The stability of 
the g-crystals is evident, since it is a triangular lattice 
that realizes the absolute minimum of the energy of mu- 
tual particle repulsion if the particles are free to move. 
With a view to formulating the general approach to the 
problem, we start with consideration of small vicinities of 
p q . As infinitesimal transformations of the system with 
the Hamiltonian ([!]) are impossible, a small change in N 
results only in formation of isolated defects in a ^-crystal, 
the space structure of the defects essentially depending 
on whether they arise on an increase or on a decrease in 
N. This fact can formally be expressed by the identity 

E g (N±SN,JV)~E g (N,M)=±p q h SN, (2) 
SN/N -> 0, 5N/N 1/2 -> oo , (3) 

where E g is the ground-state energy, SN is an arbitrary 
macroscopically small change in N. The proportionality 
coefficients, fig, are functions of defects formation ener- 
gies, interval [p~ , fig] being the range of the values of the 
chemical potential, fi, at which the g-crystal exists. 

At first glance it would seem that SN should be 
identified with a number of zero-dimensional defects 
("idefectons"), each of them arising as one particle 
is removed from g-crystal (— defecton) or added to it 
(+defecton). In such a case pf equals the energy of 
idefecton formation, ±e^. The crux of the matter is 
that this seemingly evident statement is incorrect, as 
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a rule, due to a coalescence of defectons of the same 
"sign". What this means is removing from or adding 
to a g-crystal two particles brings about a bound state of 
—defectons or +dcfectons (" bidefecton" ) whose energy is 
less than 2e^ (Fig. |l|). We have found out the coales- 
cence by computer calculations, which were carried out 
for v(r) = exp(-r/R)/r p (0 < R < 20, 1 < j3 < 12) 
and p > 10~ 2 . They have shown that the phenomenon 
occurs at any q > 2 and for all v(r) under consideration. 
The algorithm of the calculations is based on the dipole 
description of defectons which is as follows. 

Each particle displacement in a g-crystal resulting from 
idefectons formation can be considered a "dipole" of 
a sort. This is a pair consisting of the particle shifted 
and the " antiparticle" located at the host-lattice site left 
by the particle, the antiparticle "charge" being equal by 
magnitude but opposite in sign to that of the particle. A 
idefecton thus comprises a hole (— defecton) occupying 
one of the host-lattice sites, r~, or an added particle 
(+defecton) placed in the g-crystal at a free host-lattice 
site, r + , and a system of dipoles, the energy of idefecton 
formation having the form 

4 = min ^±(^) ±J2^i\n - ^1) + su^j , (4) 

where Agf(r) = f(r — £) — f(r) (f(r) is an arbi- 
trary function); f% and £j are the i-th dipole antiparticle 
site vector number and displacement vector respectively; 
±A^v(\fi — r ± |) is the energy of interaction between the 
hole /added particle and z-th dipole; n is the number 
of dipoles; ~) = — 2it g is the energy of hole forma- 
tion (u q > is the energy of the g-crystal per particle), 
v + (r + ) > is that of interaction between the (ideal) q- 
crystal and the added particle; SU > is the excitation 
energy of the g-crystal produced by the displacements of 
n particles: 

SU = + - A£ A^(|r 4 - r k \) . 

i — l i,k=l 
i^k 

Here v q {^) > is the excitation energy produced by dis- 
placement of a g-crystal particle by host-lattice vector £; 
the second term is the energy of the dipole-dipolc inter- 
action. Minimization in (^) is performed with respect to 
variables, f^, £j, n, and f + . 

To find the energy of formation of a defect with an ar- 
bitrary number of holes/added particles, m, (let their 
host-lattice sites be (a = l,...,m)) the first two 
terms in Eq. (Q) should be replaced by the correspond- 
ing sums over r^, and besides, the additional term, 
the energy of mutual holes/added particles repulsion, 
V Tep = lt2 a ^ b v(\r± - f^l), should be included, among 
the variables to be found as a result of the minimization 
being f ± . 



Using the dipole description, the nature of the coales- 
cence can be clarified. The dipoles of idefecton are al- 
ways arranged in such a way that they are all attracted to 
the hole/added particle. The bidefecton formation causes 
an extra attraction between the dipoles and holes/added 
particles, as compared with two isolated defectons. The 
energy gain exceeds y re p and thereby results in the co- 
alescence. This fact is illustrated in Fig. [|. It should 
be noted that for all q > 2 the structures of the attrac- 
tion cores of the bidefectons (the holes and their nearest 
dipole surroundings) are much the same. In the case 
q = 1 the dipoles do not arise (a —defecton is simply a 
hole), and hence, the coalescence does not occur. 

Formation of a bound state of three idefectons dimin- 
ishes the energy per particle removed/added still further, 
and so on to the extent of formation of an infinite strip of 
bound idefectons, which is a periodic ID structure with 
a primitive translation vector qa (see Fig. |l|). Hence, it 
is ID defects, but not the defectons, which are bound to 
arise in a g-crystal as a result of a small variation in N. 

The simplest ID defect in a g-crystal is a strip of rar- 
efaction or compression which arises when a g-crystal 
part adjacent to a particle line with some primitive trans- 
lation vector, d, is shifted as a whole relative to the 
other one by a host-lattice translation vector, £. For- 
mation of N s strips results in the change, ±5Af, in the 
dimensionless volume of the system: i<W = ±aLN s , 
(a = 2||* x d|/V3, L - N 1 / 2 is the length of a strip). 
If Ns/N 1 / 2 — > 0, N s — > oo j the change in the system 
energy, SE, is proportional to SAf: 

6E = e(d,Z)6Af 
s(d,0 = ^J2f'j:Zi^M\r-ne\). 1 > 

Here the accent at S means summation over the half- 
plane r = Id + me, (I = 0, il, . . . , — oo < m < 0); 
e is any of the g-crystal primitive translation vectors 
differing from d. Only such strips can be relevant to 
the ground-state properties whose vector parameters, 
d and £, realize the minimum of SE, which equals 
E g (N,Af — 5N) - E g (N,N) (rarefaction, e < 0) or 
E g (N,JV+5JV)-E g (N,J\f) (compression, e > 0). We will 
call them '—strips' or '+strips', depending on whether 
they result from rarefaction (sign — ) or compression (sign 
i); their vectors d ( 'directors') and £ are denoted here by 
dg and £^ respectively. Our computer calculations based 
on (||) have shown that the space structure of istrips is 
the same for all g and v(r): 

d q = qa a , £± = =F«/3, a ^ (3 . (6) 

where vectors a a (a = 1,2,3; \a a \ = 1) are the host- 
lattice primitive translation vectors which form an angle 
of 120° with one another. 

As indicated later, istrips of the same sign repel each 
other at any distances. Hence, the complex strip which 
result from the defecton coalescence is bound to dis- 
sociate to its constituents, simple ones. Noticing also 
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that breaking istrips into fragments shifted relative to 
one another increases inevitably the energy of the sys- 
tem (this physically evident statement is completely con- 
firmed by our computer calculations), we can conclude 
that it is just the istrips which are the above-mentioned 
ID defects. 

As follows from simple thermodynamic considerations, 
the coefficients p^ in Eq. (||) are related to the energy of 

istrip formation, ej: = \e(d^,&)\, by the expression 
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Lest there be a contradiction to the fact of the defec- 
ton coalescence, the energies of idefectons formation are 
bound to be external to the interval [p~ ,p q ]. This can 
be established by a rather simple reasoning for v(r) which 
goes to zero over distances R «f ~ q. As = 1, from 
Eq. (|5|) we have ~ (qR)~ 1 u q . On the other hand, 
|e_| is comparable with u q , and hence, p,~ ~ (q/R)u q 
|e_|. To evaluate e+ one should take into account that 
+defecton formation significantly decreases the least of 
the interparticle distances, r m i n , as compared with r. 
Therefore, e+ ~ v(r m in) ^ ~ (q/R) v (^)- To make 
sure that the inequalities, \e~\ < ji~ < pj < e q , hold 
irrespective of a form of v(r) the computer calculations 
of p q based on Eq. (Q), (||), and (Q) have been carried 
out in parallel to the computer studies of the coalescence. 
They have confirmed that the above inequalities are re- 
ally the case for all v{r) under consideration and for all 
q > 2. 

The above reasoning shows that at \p — p q \ <C 1 the 
2D LSLRIR ground-state space structure is a kind of a 
dilute solid solution of —strips or +strips (depending on 
the sign of p — p q ) in the g-crystal, istrips sharing their 
primitive translation vector, qa a , with 'unperturbed' q- 
crystal strips. It should be noted that due to the degener- 
ation caused by the hexagonal symmetry of the g-crystal 
(see Eq. (^|)) at a given p superstructures can exist which 
comprise strips with different . However, these can- 
not be the ground state, as intersections of the strips 
inevitably result in an extra increase in the system en- 
ergy. It can also be shown that the same holds for the 
mixtures of istrips with the same director but different 
displacement vectors, £ q = ^fap. To find the arrange- 
ment of the istrips with a given director and a given 
displacement vector one can consider the 2D LSLRIR as 
an ID LSLRIR whose 'particles' arc particle lines of the 
2D system, 

f=kqa a + la,3 (k = 0, il, ±2, . . . ; a ^ (3), (8) 

with integers I as 'coordinates'. From this point of view 
a istrip is a pair of lines (^) which are shifted relative to 
one another by vector (q =p l)a/5, while the relative dis- 
placement vector of the lines which form an unperturbed 
strip is qap. 

The line arrangement is governed by the pair potential 
of line- line interaction, 



V q (l 2 -h)= Yl v(\kqa a + {h-h)a \), (a^/3), 

(9) 
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where h,l2 are 'coordinates' of interacting lines. Sim- 
ple computer calculations show that due to convexity of 
v(r) the pair potential (|^) is positive and monotonically 
decreasing function of / which also meets the convexity 
condition, F q {l) = V q {l) - 2V q (l + 1) + V q (l + 2) > 0, 
over the whole range I > 1 and for all q > 2. It fol- 
lows herefrom that the arrangement of the lines (JsJ) , at 
least in some vicinity of p g , is the same as that of the 
ID LSLRIR particles- and hence, it obeys the universal 
Hubbard algorithmH~Ll 



l m = [Zm], I = ■&/ q 



(10a) 



where I is the mean line separation, = l m is the 
'coordinate' of m-th line. The algorithm determining 
the arran gement of istrips can be found directly from 
Eq. (pal): 



n m = [m/<*]; q{q T c±)=d < c± < 1. (10b) 

Here m enumerates the istrips, n m is the number of 
the particle lines between m-th istrip and istrip with 
rn = 0; = /Ni is the concentration of the istrips 
(Ng and Ni are total numbers of the istrips and the 
particle lines respectively). 

The 2D LSLRIR ground-state space structures ( |l0b| ) 
are stable not only in some vicinity of $ = q 2 , but also 
over the whole range of istrip concentrations, < < 

1. This is due to the following closely correlated facts, 
which we (using computer calculations) have found to 
take place for any v(r) under consideration: a) zero- 
dimensional defects (defectons) of a structure described 
by Eq. (|l0| ) coalesce irrespective of a value; b) there- 
fore, an infinitesimally small variation in cj produces 
only ID defects in the structure; c) their primitive trans- 
lation vectors turn out to be equal to qa a for any . 

At cj = 1 (i? = q(q+l)) algorithm ( |l0| ) determines the 
stable simple crystal with primitive translation vectors 
qa a , (q+ l)a/3 (a ^ (3). This can also be considered to be 
the structure arising at c+ = 1 as a result of compression 
of the q+ 1-crystal. Hence, universal algorithm ( [iTj| ) (q 

2. '.') .... covers the whole range 2 < d < oo, expressions 
@ holding over the intervals q(q — 1) < 1? < q(q + 1). 
For each interval there is its own set of directors, qa a . 

The only exception, as our computer calculations have 
shown, is a vicinity of d = 2 where the above algorithm 
fails in the case of slow-diminishing pair potentials v(r). 
Particularly, if v(r) = exp(r/i?)/r or v(r) = 1/r 7 , this 
occurs for R > R « 0.38 or 7 < 70 « 4.3 respectively. 
The matter is that for v(r) of such a type the coales- 
cence of idefectons of the simple crystal with -d — 2 (its 
primitive translation vectors are a ai 2afj, (a ^ (3)) leads 
to formation of ID defects the primitive translation vec- 
tor of which, 4a a , is doubled in length as compared with 
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director corresponding to q = 2. These defects cannot 
disassociate into simple strips of the above type. Our pre- 
liminary results suggest that on interval the 2 < ■& < 4 
there exists a critical point (its position depends on a 
form of v(r)), at which the 2D LSLRIR ground state un- 
dergoes a first-order structure transition with doubling 
of a spacing, 2a a — ► 4a Q . We are going to carry out a 
detailed study of the situation in a separate paper. 

The 2D LSLRIR ground-state dependence of the chem- 
ical potential, /i, on p, just as it takes place for the ID 
LSLRIRO lI, is a devil staircase whos e ste ps occur at all 
rational Taking into account Eq. (10a), one can rep- 
resent •& in a form § = ql (q = 2, 3, . . .), where I = Q/P, 
and P, Q are any coprime integers meeting conditions 
3-1 < Q/P < 9+1- Thus, the 2D LSLRIR ground state 
is a set of crystals with P particles per cell and primitive 
translation vectors qa a ,Qap, a ^ (3, the crystals with 
one particle per cell having ■& = m?. or q(q +1). Using 
results of the ID LSLRIR theorytrQ, we have found the 
widths of the devil-staircase intervals, A/i, to be related 
to the line pair potential (m by an expression 



Ap = Q mF q {Qm - 1) > 



(11) 



m— 1 



At a given q > 2 the expression holds over the interval 
q{q — 1) < $ < q(q + !)■ It should be noted that Ap, 
similarly to the ID case, does not depend on the number 
of particles per cell, P. In the case of g-crystal (Q — 
q, P = 1) expression ( |TT| ) coincides with p+ — fi~ , (see 
Eq. (0). 

The above theory holds over the interval 1 > p > 1/2 
if one applies it to the holes (empty host-lattice sites) as 
the particles, the filling factor of holes, ph, being equal 
to 1 — p. It is evident that the hole space structure cor- 
responding to a given ph is equivalent to the particle one 
with p = 1 — ph- 

In essence, the effective lowering of dimension and the 
coalescence of defectons underlying it are caused by dis- 
creteness of the system under consideration. Therefore, 
the phenomena should be expected to exist irrespective 
of the geometry of a host lattice. Our preliminary stud- 
ies have confirmed this suggestion, yet they have shown 
that for an arbitrary host lattice the structure of the 2D 
LSLRIR ground state turns out to be more sophisticated 
than in the case of a triangular host lattice. We also 
have revealed that at non-zero temperatures the lower- 
ing of dimension leads to an interplay of a strongly corre- 
lated liquid of the thermally fractured ID defects and an 
ideal gas of defectons mentioned. This results in a very 
interesting low-temperature thermodynamics of the 2D 
LSLRIR. We are going to publish the results concerning 
these issues in the near future. 

We would like to thank M. Pepper for bringing his ex- 
perimental results to our attention. 
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FIG. 1. The coalescence of —defectons in g-crystal with 
q = 4. Here o denotes host-lattice sites, • - particles; - 
antiparticles; ® - holes; — > - dipoles; the dotted lines mark 
off a single defecton, a bidefecton. The attraction core of the 
bidefecton consists of dipoles D\, Di, and holes Hi, Hi- The 
total energy of attraction Hi «-» D2 and Hi <-> D\ exceeds the 
energy of repulsion Hi <-> Hi- With an increase of number 
of removed particles the strip of bound defectons arises, the 
holes of the strip stringing out along a particle line passing 
through the sites Hi, Hi- 
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